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Abstract 

We have tested the theoretical values of critical exponents, predicted for 
the three-dimensional Heisenberg model, based on the published Monte Carlo 
(MC) simulation data for the susceptibility. Two different sets of the criti- 
cal exponents have been considered - one provided by the usual (perturbative) 
renormalization group (RG) theory, and another predicted by grouping of Feyn- 
man diagrams in (p 4 model (our theory). The test consists of two steps. First 
we determine the critical coupling by fitting the MC data to the theoretical 
expression, including both confluent and analytical corrections to scaling, the 
values of critical exponents being taken from theory. Then we use the obtained 
value of critical coupling to test the agreement between theory and MC data at 
criticality. As a result, we have found that predictions of our theory (7 = 19/14, 
T] = 1/10, w = 3/5) are consistent, whereas those of the perturbative RG the- 
ory (7 ~ 1.3895, 77 ~ 0.0355, ui ~ 0.782) are inconsistent with the MC data. 
The seemable agreement between the RG prediction for r\ and MC results at 
criticality, reported in literature, appears due to slightly overestimated value of 
the critical coupling. Estimation of critical exponents of 3D Ising model from 
complex zeroth of the partition function is discussed. A refined analysis yields 
the best estimate 1/v ~ 1.518. We conclude that the recent MC data can be 
completely explained within our theory (providing 1/v — 1.5 and u> = 0.5) 
rather than within the conventional RG theory. 

Keywords: Heisenberg model, Ising model, Monte Carlo data, critical expo- 
nents, partition function zeroth 

1 Introduction 

In our previous work jlj], we have reported the possible values of exact critical expo- 
nents for the Ginzburg-Landau phase transition model predicted from a reorganized 
perturbation theory. These predictions are in exact agreement with the known ex- 
act and rigorous results in two dimensions ||], and are equally valid also in three 
dimensions. Our predictions have been compared to some original data of Monte 
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Carlo (MC) simulations and experiments ||, ||, ra, and a remarkable agreement has 
been found. 

However, there is still rather paradoxical and unclear situation regarding the 
MC results. On the one hand, we have shown theoretically Q the invalidity of the 
conventional RG expansions || fj], ||], but, on the other hand, the published papers 
on MC simulations usually claim to confirm the values of critical exponents coming 
from these expansions and being in contradiction to our results. 

Contrary to the usual claims in the published papers that the values of critical 
exponents can be obtained from the Monte Carlo data with a striking accuracy, i. e. 
with an error much smaller than 0.01, our expierence in analysis of several such data 
shows that in reality it is very difficult to extract accurate and reliable estimates 
therefrom. The problem is that a fitting of MC data to a simple theoretical ansatz 
(including no corrections to scaling) can provide a rather small statistical error, but 
the obtained result is not reliable since it contains an uncontrolled systematical error 
due to the neglected corrections to scaling. Moreover, confluent (i. e., those related 
to the universal properties of the critical point) and analytical corrections can be 
equally important at finite values of the reduced temperature at which the simula- 
tions have been done, since the amplitude of the leading analytical correction can 
be remarkably larger than that of the confluent correction. Our analysis of the sus- 
ceptibility data for the three-dimensional Heisenberg model (Sec. ||) has shown that 
the estimated value of the critical exponent 7 decreases by several percents due to 
the confluent correction, and the result can be changed remarkably by the analytical 
correction too. Thus, both kind of corrections should be taken into account, but 
this is not possible in the usual applications related to the determination of critical 
exponents, since inclusion of both kind of corrections in a theoretical ansatz strongly 
increases the statistical errors. 

As regards the fitting of MC data at criticality, only confluent corrections are 
present, but the usual estimations are rather sensitive to the precise value of the 
critical coupling. In this aspect, our estimation of the critical exponent rj Q from 
the MC simulated fractal dimensionality of the three-dimensional Ising model at the 
critical point (i. e., from MC data of Q) is preferable to a more conventional, but 
much more sensitive to the precise value of the critical coupling f3 c , estimation of 
this exponent from the susceptibility data at criticality. According to the published 
results P], the second method seems to give smaller values of r\ in three dimensions 
(about 0.027 for Heisenberg model |J) as compared to the first one (about 1/8 for 
Ising model Q), but the reason for the discrepancy could be an inaccuracy in the 
estimated value of (3 C . In the case of the Heisenberg model, this value has been 
overestimated, indeed, as discussed in Sec. § More recent MC results reported 
in (T^j also provide rather small values of r\ (about 0.04 for 0(n) models with n = 
2, 3, 4). However, the infinite volume extrapolation in ju| is erroneous in view of our 
theory, and a selfconsistent treatment, based on our theoretical predictions, reveals 
no contradiction to the MC data (Sec. ||). 

In the present work we have proposed a Monte Carlo test, based on a high 
quality susceptibility data |J , where the above discussed problems with corrections 
to scaling are solved on a higher level than it has been done in the currently published 
papers. Namely, our method enables us to test the agreement of MC data with given 
(fixed) theoretical values of critical exponents by taking into account both the leading 



2 



confluent and the leading analytical correction. Our test consists of a very accurate 
determination of the critical coupling followed by a fitting of the susceptibility data 
at criticality. It has shown (Sec. ||) that the actually discussed MC data are in 
agreement with our theoretical values of critical exponents, but not with those of 
the RG expansions. 



2 Critical exponents from our theory 

Our theory provides possible values of exact critical exponents 7 and v for the tp 
model whith 0(n) symmetry (n-component vector model) with the Hamiltonian 

H/T = J [?w 2 (x) + c(V(/?(x)) 2 + V(x)] dx , (1) 

where tq is the only parameter depending on temperature T, and the dependence is 
linear. At the spatial dimensionality d = 2, 3 and n = 1, 2, 3, ... these values are [[jj 

d + 2j + 4m 

7 d(l + m + j) -2j ' U 
_ _ 2(1 +m)+j 



d(l + m + j)- 2j 



where m > 1 and j > —m are integers. At n = 1 we have m = 3 and j = to fit the 
known exact results for the two-dimensional Ising model. As proposed in Ref. [jij, 
in the case of n = 2 we have m = 3 and j = 1, which yields in three dimensions 
^ = 9/13 and 7 = 17/13. 

In the present analysis the correction-to-scaling exponent 9 for the susceptibility 
also is relevant. The susceptibility is related to the correlation function in the Fourier 
representation G(k), i. e., x °c G(0) [f7|. In the thermodynamic limit, this relation 
makes sense at T > T c , where T c is the critical temperature. According to our 
theory, G(0) can be expanded in a Taylor series of t 2v ~^ at t — > 0. In this case the 
reduced temperature t is defined as t = tq(T) — ro(T c ) oc T — T c . Formally, t 2 ' y ~ dv 
appears as second expansion parameter in the derivations in Ref. |IJ, but, according 
to the final result represented by Eqs. @) and (H|), (27 — du)/(2u — 7) is a natural 
number. Some of the expansion coefficients can be zero, so that in general we have 

9 = l(2v-j), (4) 

where t may have integer values 1, 2, 3, etc. One can expect that t = 4 holds at 
n = 1 (which yields 9 = 1 at d = 2 and 9 = 1/3 at d = 3) and the only nonva- 
nishing corrections are those of the order t e , t 2e , t 3e , since the known corrections 
to scaling for physical quantities, such as magnetization or correlation length, are 
analytical in the case of the two-dimensional Ising model. Here we suppose that 
the confluent corrections become analytical, i. e. 9 takes the value 1, at d = 2. 
Besides, similar corrections to scaling are expected for susceptibility \ and mag- 
netization M since both these quantities are related to G(0), i. e., % oc G(0) and 
M 2 = lim a; _ >oo ((/3(0)<^(x)) = limy^oo G(0)/V hold where V = L d is the volume and 
L is the linear size of the system. The above limit is meaningful at L — > 00, but 
G(0)/V may be considered as a definition of M 2 for finite systems too. The latter 
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means that corrections to finite-size scaling for x an d M are similar at T = T c . 
According to the scaling hypothesis and finite-size scaling theory (Sec. ||), the same 
is true for the discussed here corrections at t — * 0. Thus, the expected expansion of 
the susceptibility x looks like x = ( a o + CL±t e + d2t 2e + ■ ■ -j . 

Our general hypothesis is that j = j(n) and I = £(n) monotoneously increase 
with n to fit the known exponents for the spherical model at n — > oo. In particular, 
we expect that j(n) = n — 1, £(n) = n + 3, and m = 3 hold at n = 1,2,3 and, 
probably, also in general. This hypothesis is well confirmed by MC results discussed 
here and in Ref . H . 

We allow that different t values correspond to the leading correction-to-scaling 
exponent for different quantities related to G(k). The expansion of G(k) by itself 
contains a nonvanishing term of order i 2l/ ~ 7 = t* 1 " (in the form G(k) ~ 
t" 7 [g(kt~ u ) + t rru g 1 (kt~ v )] whith g^O) = 0, since i > 1 holds in the case of suscep- 
tibility) to compensate the corresponding correction term (produced by c(V</?) 2 ) in 
the equation for 1/G(k) (cf. Q). The latter means, e. g., that the correlation length 
£ estimated from an approximate ansatz like G(k) oc 1/ [k 2 + (l/£) 2 ] used in 1C] 
also contains a correction proportional to W . Since r]v has a rather small value, 
the presence of such a correction (and, presumably, also the higher order correc- 
tions t 2l]U , t 3r,L/ , etc.) makes the above ansatz unsuitable for an accurate numerical 
correction-to-scaling analysis. Due to this reason the susceptibility data, but not 
the correlation length data of Ref. Q, are used in our further analysis. 

The correction t 71 " is related to the correction L _r? in the finite-size scaling expres- 
sions at criticality. The existence of such a correction in the asymptotic expansion of 
the critical real-space Green's (correlation) function, i. e. 
G(rL) oc L 2 ~ v ~ d (1 + o(L~ v )) where r is a constant, is well confirmed by our recent 
(preliminary) results for the 2D Ising model. These results for L = 2, 4, 6, . . . , 16 
have been obtained by an exact numerical transfer-matrix algorithm. In such a 
way, the infinite volume extrapolation in [1C] appears to be incorrect, therefore the 
obtained there results do not represent a serious argument against our theory. More- 
over, if the extrapolation in |]l0| is done including the correction then the results 
for 0{n) models with n = 2,3 appear to be in a satisfactory agreement (within the 
extrapolation errors) with our values r] = 1/9 and rj = 1/10, respectively. 

Our consideration can be generalized easily to the case where the Hamiltonian 
parameter tq is a nonlinear analytical function of T. Nothing is changed in the 
above expansions if the reduced temperature t, as before, is defined by t = ro(T) — 
ro(T c ). However, analytical corrections to scaling appear (and also corrections like 
(T — T c ) m+nd with integer m and n) if t is reexpanded in terms of T — T c at T > T c , 
The solution at the critical point remains unchanged, since the phase transition 
occurs at the same (critical) value of tq. 



3 Estimation of the critical exponent 7 from MC data 

In this section we discuss the estimation of the susceptibility exponent 7 for the 
classical three-dimensional Heisenberg model. Our analysis is based on the fit- 
ting of the susceptibility (MC) data to a theoretical ansatz. According to the 
finite-size scaling theory, the susceptibility x depending on the reduced temperature 
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t = 1 — (3/ f3 c (where t > 0) and the linear size of the system L reads 

X = L^g(L/t), (5) 

where g(L/£) is the scaling function and £ ~ t~ v is the correlation length of an 
infinite system. Eq. (|5|) holds precisely at L — > oo and t — > for any given value of 
L/£. At finite values of i and L corrections to (||) exist. Eq. (||) can be rewritten as 

x = rv (tL 1 /-) , ( 6 ) 

where g(y) = y^l v f (y 1 ^^ ■ In the thermodynamic limit L —> oo Eq. (JgJ) reduces to 
X = bot -1 , where bo = lim fix) is the amplitude. A natural extension of Eq. <M), 

x—*oo 

including corrections to scaling, is 

X = t^J2 t7l fi( tLl/U ) > (7) 

where 70 = 0, /o(^) = /(^), and the terms with Z > represent all the corrections 
in the asymptotic expansion of % at i — > for any given value of x = tL l l u . In the 
thermodynamic limit we have lim fi(x) = bi, where b\ are the amplitudes. The most 

x— >oo 

important correction terms in the sum over I are the leading confluent correction 
b\ £ 71 with the exponent 71 = 9 and the leading analytical correction 62 f 12 with 
72 = 1. Although < 1 holds, the analytical correction also should be included 
at finite values of t used in practical simulations: because of absence of a direct 
correlation between the amplitudes of confluent and analytical corrections, the ratio 
r = b<ijb\ can be arbitrarily large. One can expect that the higher order confluent 
corrections (i. e., those proportional to t 2e , t 3e , etc.) are small as compared to the 
leading confluent correction, and the same is true for analytical corrections. We 
consider the case of small t and large x, i. e., small fi(x) — b[. In this case Eq. (]?]) 
can be written as 

x ~ t -y\l + b (t e + rt)] f(tL x ' v ) , (8) 



where b = bi/bo is a constant. 

We have used the susceptibility data simulated by an improved (cluster) MC 
algorithm reported in Ref. ]|] (x tmp vs (3, tab. IV in ||) to estimate the critical 
exponent 7 by fitting the data to Such an estimation has been done in Q, ne- 
glecting either the analytical or the confluent correction and setting / (tL l / v ^j = b$. 

Since in the actual simulations the scaling argument x = tL l / u has large enough 
values, about 6 or 7, which are varied only slightly, the latter approximation is rea- 
sonable. We have used even better approximation where ln/(x) has been linearized 
within the narrow range of x variation, and the simulated data points for In x have 
been fitted to the resulting theoretical expression 

]nx(t,L) = a-7lnt + ln l + b(t 6 + rt\ +ptL 1/u , (9) 

where a and p are constants. The minimum of the sum of the squared deviations for 
N data points S(N) corresponds to the least-squares fit. Besides, it is reasonable 
to use the least-squares method just for ln%, but not for x, since the errors for ln X 
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Figure 1: Estimation of the critical exponent 7 in 3D Heisenberg model. Solid line shows 
the standard deviation a of the simulated data points from the analytical curve (^) as a 
function of 7 with parameters a, &, p, and [3 C obtained from the least-squares fit at v = 5/7, 
= 3/7 (our theoretical values), and r — 0. The dotted curve corresponds to fixed p = 0. 
The minimum of the solid curve gives the least-squres estimate 7 = 1.345. All fits (for 
different data sets) lie in the marked area which is shifted only slightly, as indicated by thin 
vertical dashed lines, if the RG values of v and 8 are used instead of ours. Our theoretical 
value 7 = 19/14 (thick vertical dashed line) is inside of the marked region, whereas that of 
the RG theory (vertical dot-dot-dashed line at 7 = 1.3895) is outside. 



data points are comparable, i. e., the relative but not the absolute errors are more 
or less equal. At large N, the inaccuracy in the fitted curve due to the statistical 
errors can be characterised by the standard deviation a = (S(N)/[N(N — l)]) 1 / 2 . 
Obviously, the minimum of a corresponds to the least-squares fit at any given N. 

We have illustrated in Fig. [I] the estimation of 7 by minimizing a with respect 
to the parameters a, b, p, and (3 C (where /3 C is incorporated in (||) via t = 1 — /3/j3 c ) 
at fixed exponents 6 = 3/7 and v = 5/7, taken from our theory (Sec. ^). The 
analytical correction to scaling has been neglected by setting r = 0. The solid line 
shows the accuracy of the fit, i. e. the value of a, depending on the choice of the 
exponent 7. The minimum of a, indicated by a vertical dotted line, is located at 
7 ~ 1.345, which corresponds to the least-squares fit. The dotted curve corresponds 
to the case of fixed p = 0. From this we can see that inclusion of the term ptL 1 / 11 
in (0), responsible for the variation of the scaling function f(x), affects the result 
only slightly. 

In spite of the very high accuracy of the fit (about 0.02% error in x), the minimum 
in a is too broad for a reliable estimation of 7 with, e. g., ±0.01 accuracy. This is a 
problem which usually appears if we use a high-level approximation including many 
fitting parameters. If the analytical correction also is included, then the situation 
becomes even worse, i. e., the a vs 7 plot is an almost horizontal line. Neglection of 
both (confluent and analytical) corrections, as it has been done finally in Q, is not 
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a solution of the problem since the result is affected significantly by the confluent 
correction. Namely, the obtained value of 7 is shifted from 1.389 to 1.345. According 
to our estimation, the statistical error for the latter result is remarkably smaller than 
the difference between these two values, so that the second value is better. Another 
problem is that the estimated value of 7 depends on 9 and v. This effect, however, 



is relatively small. By the conventional RG values 9 = 0.55 and v = 0.7073 [11] we 
obtain 7 ~ 1.354. 

Like in Ref. Q, we have estimated the possible statistical error of our result 
7 ~ 1.345 by comparing the values of 7 for a large number of different data sets 
generated from the original one (with 18 data points) by omitting some (1 to 6) 
data points. The data points have been omitted more or less randomly, but not the 
neighbouring points and not the first and the last point simultaneously, to ensure 
a sufficiently uniform distribution of the used t values and to avoid a significant 
narrowing of the total interval covered by these values. The largest deviations from 
the central 7 value 1.345 have been observed omitting the data points No. 1, 6, 10, 
14, and 17 (tab. IV in |J), which yielded 7 ~ 1.322, and the data points No. 2, 5, 
8, 11, and 14, which yielded 7 ~ 1.366. Thus, all the fits gave 1.322 < 7 < 1.366 
at r = 0, 9 = 3/7 and v = 5/7. This interval is marked in Fig. |l] by thin solid 
lines. At 9 = 0.55 and v = 0.7073 the borders of this region are shifted slightly, as 
indicated by thin vertical dashed lines. These manipulations enable us to estimate 
the possible statistical error in both cases, i. e., 7 = 1.345 ± 0.023 at 9 = 3/7 and 
7 = 1.354 ± 0.020 at 9 = 0.55. These, in fact, are maximal errors, i. e., since we 
never have observed larger deviations, the probability that the value extracted from 
exact data would be outside of the error bars is vanishingly small. 

It is a remarkable fact that our theoretical value 7 = 19/14 ~ 1.35714 (thick 
vertical dashed line) lies inside the region of maximal statistical errors, whereas that 
of the RG theory, i. e. 7 ~ 1.3895 indicated by a do-dot-dashed line, is clearly 
outside of this region. This result can be changed by the analytical correction. 
However, if the ratio of amplitudes r in Eq. (^) is positive, then the least-squares 
fit with respect to the parameters a, b, p, and f3 c always yields the central value of 
7 (with all 18 data points included) in the range from 1.345 to 1.369 at 9 = 3/7 
and v = 5/7. Here 7 = 1.369 ± 0.013 corresponds to the case of purely analytical 
correction to scaling obtained by formally setting 9 = 1. In such a way, selfconsistent 
estimations at r > yield 7 values which are reasonably close to our prediction 
7 = 19/14. Precise agreement is reached at r ~ 1.17. 

Unfortunately, we have no proof that r is positive. If we allow that r < 0, then 
a large uncertainty appears. In this case 7 can take the values as small as, e. g., 1.1 
(at 9 = 3/7, v = 5/7, and r -1.35) and as large as 1.4 (at 9 = 0.55, v = 0.7073, 
and r ~ —1.8). The actual MC data do not allow to find the true value of r (unless 
we assume that our exponents are true and, therefore, r m 1.17), since the standard 
deviation of the least-squares fit is almost independent on r. 

4 Estimation of the critical coupling 

Based on the method developed in Sec. ||, here we determine the critical coupling 
j3 c for the three-dimensional Heisenberg model assuming that critical exponents 7, 
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Figure 2: Estimation of the critical coupling (3 C by fitting the MC data to ansatz ^ with 
fixed exponents taken from our (thick solid line) and RG (thin solid line) theory. The 
minimums of a vs /3 C , used as a fitting parameter, give the least-squares estimates for the 
true f3 c value, as indicated by vertical dotted lines and arrows. The vertical dashed line 
indicates a value of f3 c proposed in || [Io|. 



9, and v are known from theory. The latter ensures a very small statistical error. 
The coefficients in (^) are found by the least-squares method. The resulting value 
of the standard deviation a vs /3 C , used as a fitting parameter, is shown in Fig. ^. 
The thick solid line corresponds to our critical exponents 7 = 19/14, v = 5/7, 
and 9 = 3/7, whereas the thin solid line - to the conventional (RG) exponents 
7 = 1.3895, v = 0.7073, and 9 = 0.55. The minimums of these curves, indicated 
by vertical dotted lines, are located at (3 C ~ 0.692795 and (5 C ~ 0.692855, respec- 
tively, corresponding to the least-squares estimates for the true values of the critical 
coupling. The estimation of maximal statistical errors, as in Sec. ||, leads to the 
following conclusions: 

1. If our values of the critical exponents 7 = 19/14, v = 5/7, and 9 = 3/7 are 
true, then 

& = 0.692795 +a00003 °. (10) 
-0.000043 

2. If the true values of the critical exponents are close to those predicted by the 
RG theory, i. e., 7 = 1.3895, v = 0.7073, and 6 = 0.55, then 

+0.000029 . , 

Pc = 0.692855 . (11) 

-0.000043 

The estimation in || gave j3 c ~ 0.6930. This value is indicated in Fig. || by 
a vertical dashed line. As we see, it clearly does not correspond to the best fit. 
To clear up the reason for the discrepancy, let us discuss the Binder's cumulant 
crossing technique used in Q and [10] for the estimation of /3 C . In this approach, the 
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Figure 3: Estimation of the critical coupling by the Binder's cumulant crossing technique. 
The straight line represents the least-squares fit of ([l2]) to the MC data for crossing points 
T cross . The zero intercept gives /3 C = 1/T C ~ 0.69286. For comparison, the approximation 
Tcross —T c oc l/ln& (where b is Binder parameter) is shown by thin solid line. 

magnetization cumulants for different lattice sizes L and V are plotted as a function 
of (5 to find the intersection point = f3 cross . According to the theory ]l2| , 

(3 cross (L, b) - p c oc L-d/")-" _ 7 __ (12) 

holds at large L, where L is the size of the smaller lattice, b = L'/L is the Binder 
parameter, and u = Qjv. The estimation in ||] has been done by approximating 
the term p = (1 — b~ w ) / (b 1 ^ — 1) in Eq. ([12] ) with const/ In b. We have made and 
have illustrated in Fig. |3] our own estimation of C from Eq. (|l2"|), using the data 
for T cross = l/Pcross extracted from Fig. 3 in Ref. ||. According to (|P2|), f3 cr0 ss is 
a linear function of p at a fixed L. The same is true for T cross in vicinity of T c . 
The straight line in Fig. |3| corresponds to the linear least-squares fit for T cross vs 
p at L = 16 (with our exponents uj = 3/5 and v = 5/7) which yields (at p = 0) 
P c ~ 0.69286. This value agree with @ and (|TTJ) within the error bars ±0.0001 
proposed in Q. In fact, the data points are too much scattered to consider such an 
estimation reliable. Due to this reason, we have not tried to estimate (3 C from the 
data of L = 12 which are even more scattered. We have depicted in Fig. |3] by thin 
solid line the fit, corresponding to the approximation T cross — T c oc 1/ In b, made in M 
at L = 16. As we see, in the scale where the original ansatz ([E|) yields a straight 
line this approximation is represented by a curve providing an underestimated value 
of T cross at p = 0, i. e., an overestimated j3 c ~ 0.6930 instead of (3 C ~ 0.69286. 
Obviously, this approximation is the reason for the discrepancy. Note that other 
kind of estimations in || provided a bit smaller (3 C values, closer to ours. 
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Figure 4: Our fit to the susceptibility (x) data of [^) for 3D Heisenberg model at and near 
criticality. Only two coefficients a and b have been used as fitting parameters in (|ll|) for 
the thick solid curve and solid circles representing In (x/L 2 ) vs InL at criticality according 
to our critical exponents (77 = 0.1, 7 = 19/14, uj = 0.6) and /3 C ~ 0.692795 estimated 
independently in Sec. |J. The same fit at RG exponents with the corresponding (3 C value 
0.692855 is represented by the dot-dot-dashed line and empty circles. Thin solid lines show 
our three-parameter fit at (3 = 0.6925 (empty squares), /3 = 0.6930 (crosses), and (3 — 0.6933 
(empty rhombs). The linear fit of |5) is shown by tiny dashed line. 



5 The test of consistency at T = T c 

Consequently following the conclusions flip] ) and (11) made in Sec. here we test 
the agreement between theory and MC data at criticality. 

According to the finite-size scaling theory, the susceptibility at the critical point 
is given by 

X <xL-r/ v (l + bL- u + ...) , (13) 

where b and u = 9/v are the amplitude and the exponent of the leading correction 
to scaling. The dots stand for further corrections. Eq. (|l^) can be rewritten as 



In 



(x/L 2 ) ~ a - 7] In L + In (l + bL~ w ) , (14) 



where a is a constant and r] = 2 — 7/1/ is the critical exponent describing the 
asymptotic long-wave behavior of the correlation function (i. e. G(k) ~ k~ 2+T> ) at 
T = T c . We have read from Fig. 6 in Ref. |J the values of \ near j3 c and have 
made the linear interpolation between j3 = 0.6927 and (3 = 0.6929 to estimate x a t 
the values of the critical coupling given by ([K]) and (pT|). So obtained x values are 
depicted in Fig. ||| by solid and empty circles, respectively. The corresponding two 
parameter (a and b in Eq. (|l4"|)) least-squares fits with fixed exponents are shown 
by thick solid line (our case) and dot-dot-dashed line (RG case). If r\ is considered 
as a fitting parameter, then in our case the least-squares fit yields rj ~ 0.105 in close 
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agreement with the theoretical value 0.1, whereas in the RG case it yields rj ~ 0.076 
in a remarkable disagreement with the theoretical value 0.0355. It can be seen 
also from Fig. || that the dot-dot-dashed line with a = —0.17034 and b = 0.1178, 
obtained at fixed «t) = 0.0355, does not provide a satisfactory fit to the data, i.e., 
this line is curved in a wrong direction. 

Our values of critical exponents provide an excellent fit to the MC data not 
only at (3 = /3 C , but also at small deviations t = 1 — j3/(3 c from the critical point 
considered in Fig. 6 of Ref. g. Our fit x = 1.1266 L 1 - 9 (1 - 0.4944L" - 6 - 0.58 tL 1A ) 
is shown in Fig. ^| by solid lines. This approximation is consistent with the finite-size 
scaling theory at large L and small tL 1 ^. The data points in Fig. |3] correspond to 
12 < L < 48. At smaller L values the second-order corrections to scaling, neglected 
in our ansatz, could be relevant. The solid curve at /3 = 0.6930 is the most linear 
one within 12 < L < 48, as it is evident from Fig. [| where the straight -line fit of || 
is shown by a tiny dashed line. It is evident also that the good linearity of In (x/£ 2 ) 
vs InL in this region does not mean that f3 c ~ 0.6930 and r/ ~ 0.027. 

One of the arguments in || , supporting the idea that 77 has a very small value 
(rj < 0.05), is based on the simulated data for \ vs the correlation length £ for finite 
systems. However, the variation of r\ e ff with £ in Fig. 12 of Ref. H can be well ex- 
plained by presence of corrections of the kind £~ W7? , where m = 1, 2, ... and rj = 1/10, 
consistent with the correction-to-scaling analysis in Sec. |2| (see remarks regarding 
the actual approximation for £). In this case rj e ff can behave nonmonotoneously, as 
well. As regards other arguments in || in support of the conventional RG values of 
critical exponents, they are weaker than our contraarguments discussed here, since 
all the final estimates in Q are obtained neglecting corrections to scaling. Note also 
that such kind of simple estimations not always give very small values of r/. In par- 
ticular, the values of about 0.15 follow from MC study of Heisenberg fluid O]. In 
view of our theory and presented here analysis of the MC method, the discrepancy 
between the so called "lattice" and "off-lattice" critical exponents discussed in |l3j 
can be well understood as an error of about ±0.07 (in 77) of the above discussed 
simple estimations. 

It is noteworthy that a large variety of experimental measurements in Ni dis- 
cussed in [|l4|], confirm our values of critical exponents 7 = 19/14 = 1.357..., j3 = 
(d-2+n) u/2 = 11/28 = 0.3928..., and 5 = (d+2-ry)/(d-2+r/) = 49/11 = 4.4545... 
rather than those of the RG theory (7 ~ 1.3895, 5 ~ 4.794, and j3 ~ 0.3662 (nfl). 



6 Comparison to 3D Ising model 



In this section we discuss the recent MC results [15] for the complex zeroth of the 
partition function of the three-dimensional Ising model. Namely, if the coupling (3 
is a complex number, then the statistical sum has zeroth at certain complex values 
of j3 or u = e~@. The nearest to the real positive axis values /3j and u® are of special 
interest. Neglecting the second-order corrections, u® behaves like 

u? = u c + AL- l > u + BL- { - l l^- ul (15) 

at large L, where u c = e~@ c is the critical value of u, A and B are complex constants, 
and uj is the correction-to-scaling exponent. According to the known results (see, 
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e. g., the solution given in |]16|), the partition function zeroth correspond to complex 
values of sinh(2/3) located on a unit circle in the case of 2D Ising model, so that A is 
purely imaginar. The latter means that the critical behavior of real and imaginary 
parts of u\ — u c essentially differ from each other, i. e., Re (u\ — u c ) oc L~( 1 / !y )~ a -' and 
Im (ui) oc L~ x l v (where, in this case of d = 2, v = uj = 1) at L — > oo. The MC data 
of [15 1, in fact, provide a good evidence that the same is true in three dimensions. 

Unfortunately, the authors of Ref. [[D| have not tried to find the objective truth 
regarding the behavior of the complex zeroth, but only have searched the way how 
to confirm the already known estimates for v. Their treatment, however, is rather 
doubtful. First, let us mention that, in contradiction to the definition in the paper, 
u\ values listed in Tab. I of [[D| are not equal to e - ^ 1 (they look like e -4 ^ 1 ). Second, 
the fit to a theoretical ansatz for | u®(L) — u c |, Eq. (6) in fl5| , is unsatisfactory. This 
ansatz contains a misterious parameter 03. If we compare Eqs. (5) and (6) in [15], 
then we see immediately that 03 = (l/v) + uj. At the same time, the obtained 
estimate for 03, i. e. 03 = 4.861(84), is completely inconsistent with the values of 
(l/^)+a>, about 2.34, which follow from authors own considerations. Our prediction, 
consistent with the correction-to-scaling analysis in Sec. [2| (and with t = 4 in (|2|) 
to coincide with the known exact result at d = 2), is v = 2/3 and uj = 1/2, i. e., 
(l/i/) + w = 2. 

To obtain a more complete picture, we have considered separately the real part 
and the imaginary part of — u c . We have calculated from (3® data listed 
in Tab. I of Jl^] and have estimated the effective critical exponents y' e ^(V) and 
y"fj(L), separately for Re (u® — u c ) and Im (u^), by fitting these quantities to an 

ansatz const ■ L~ Ve ff and const ■ L~ Ve ff, respectively, at sizes L and L/2. The 
value of u c consitent with high- and low- temperature series as well as MC |l8| 
estimations of the critical coupling, f3 c ~ 0.221659, have been used. The results 
are shown in Fig. ||. As we see, y' e ^ (empty circles) claims to increase above y" e ^ 
(solid circles) when L increases. This is a good numerical evidence that, like in the 
two-dimensional case, the asymptotic values are y' = lim^oo y' e ^{V) = (1/^) + uj 
and y" = lim.L-,ooy"ff(L) = 1/v. According to our theory, the actual plots in the 

L -1 / 2 scale are linear at L — > 00, as consistent with the expansion in terms of . 
The linear least-squares fits are shown by solid lines. The zero intercepts 1.552 and 
1.913 are in approximate agreement with our theoretical values 1.5 and 2 indicated 
by horizontal dashed lines. The relative discrepancy of about 4%, presumably, is 
due to the extrapolation errors and inaccuracy in the simulated data. 

The behavior of y' e ff is rather inconsistent with the RG predictions. On the one 
hand, y' e ^ claims to increase above y" e ^ and also well above the RG value of \jv 
(the lower dod-dot-dashed line at 1.5863), and, on the other hand, the extrapolation 
yields y' value (1.913) which is remarkably smaller than (l/v) + uj ~ 2.3853 (the 
upper dot-dot-dashed line) predicted by the RG theory. For selfconsistency, we 
should use the linear extrapolation in the scale of L~ u with uj = 0.799 (the RG 
value). However, this extrapolation (tiny dashed line in Fig. yielding y' ~ 1.757, 
does not solve the problem in favour of the RG theory. 

The data points of y' e ^ look (and are expected to be) less accurate than those of 
y'eff 1 since Re («? — u c ) has a very small value. The y" e ^ data do not look scattered, 
therefore they allow a refined analysis with account for nonlinear corrections. To 
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Figure 5: Effective critical exponents for the real (empty circles) and the imaginary (solid 
circles) part of complex partition-function-zeroth of 3D Ising model depending on L -1 / 2 , 
where L is the linear size of the system. Solid lines show the linear least-squares fits. The 
asymptotic values from our theory are indicated by horizontal dashed lines, whereas those 
of the RG theory - by dot-dot-dashed lines. A selfconsistent extrapolation within the RG 
theory corresponds to the tiny dashed line. 



obtain stable results, we have included the data for smaller lattice sizes L = 3 and 



L = 4 given in [19]. In principle, we can use rather arbitrary analytical function 



4>((3) to evaluate the effective critical exponent 

y'i ff {L) = In [lm<t> (#(L/2)) /Im (ft (#(£))] / In 2 

and estimate its asymptotic value y" . For an optimal choice, however, y" e ^{V) vs 
L~ u plot should be as far as possible linear to minimize the extrapolation error. In 
this aspect, our choice (ft = exp(— 0) is preferable to (ft = exp(— 4/3) used in fl9[| . We 
have tested also another possibility, i. e. (ft = sinh(2/3), providing almost optimal re- 
sults in the case of 2D Ising model. In Fig. ^ we have shown the slope of y'lff vs L -1 / 2 
curve, calculated from the MC data of [[H| [HJ, for (ft = exp(— /?) (empty circles) and 
(ft = sinh(2/3) (solid circles). It is evident that in both cases the slope cannot be rea- 
sonably approximated by a linear function of L -1 / 2 , but can be quite well described 
by a parabola. The latter means that y" e ^(V) can be satisfactory well approximated 
by a third-order (but not by a second-order) polinomial in L -1 / 2 . The correspond- 
ing four parameter least-squares fits are shown in Fig. 0. They yield y" ~ 1.473 
in the case of (ft = exp(— j3) (long-dashed line) and y" ~ 1.518 at (ft = sinh(2/3) 
(solid line). It is evident from Fig. |7] that in the latter case we have slightly better 
linearity of the fit, therefore 1/u ~ 1.518 is our best estimate of the critical exponent 
1/u from the actual MC data. Thus, while the row estimation provided the value 
y" = 1/f ~ 1.552 which is closer to the RG prediction 1/v ~ 1.5863 (horizontal 
dot-dot-dashed line), the refined analysis reveals remarkably better agreement with 
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Figure 6: Slope of the y"ff vs L x / 2 plot in Fig. || (including also smaller sizes L). The 
empty circles correspond to <j> = exp(— /3), whereas the solid circles to <j) = sinh(2/3). The 
corresponding least-squares fits 1.1840 - 8.507L~ 1/2 + f and 0.6669 - 3.645L" 1 / 2 + 

9.275L" 1 are shown by long-dashed line and solid line, respectively. 




Figure 7: Effective critical exponent y'^j(L) for the imaginary part of the complex 
partition-function-zeroth as a function of L -1 / 2 , where L is the linear size of the system. 
The empty circles correspond to 4> = exp(— f3), whereas the solid circles to <f> — sinh(2/3). The 
corresponding least-squares fits 2/"//( L ) = L4731 + 1-3345L" 1 / 2 - 4.7657X- 1 + 6.8962£- 3 / 2 
and y" ff (L) = 1.5180 + 0.7301iT 1/2 - 2.0397L" 1 + 3.3181i" 3/2 are shown by long-dashed 
line and solid line, respectively. Our asymptotic value y" = 1/v = 1.5 is indicated by 
horizontal dashed line, whereas that of the RG theory (1.5863) - by dot-dot-dashed line. 
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our (exact) value \ jv = 1.5 (horizontal dashed line). 

7 Conclusions 

In summary, we conclude the following. 

1. Corrections to scaling for different physical quantities near and at criticality 
have been discussed in framework of our recently developed theory |Q] (Sec. |2|). 

2. The critical exponent 7 for 3D Heisenberg model has been estimated by fitting 
the original susceptibility (MC) data of || to an ansatz of finite-size-scaling 
theory which includes the leading confluent correction-to-scaling term (Sec. |3|). 
The obtained estimates (7 = 1.345±0.023 and 7 = 1.354±0.020) agree within 
error bars with our theoretical value 7 = 19/14 ~ 1.35714 and disagree with 
the conventional RG value 7 ~ 1.3895. Taking into account also the leading 
analytical correction, a selfconsistent estimation always yields the central value 
of 7 in the range of 1.345 < 7 < 1.369, i. e., reasonably close to our (exact) 
value 19/14, if the ratio of amplitudes r for analytical and confluent corrections 
is varied from to 00. 

3. Based on MC data for susceptibility in 3D Heisenberg model, a very accurate 
estimation of the critical coupling has been made at given values of critical 
exponents (Sec. |I|), taking into account both confluent and analytical correc- 
tions to scaling. These estimates, combined with fits in vicinity of the critical 
point (Sec. |5|), allowed us to test the consistency between theoretical values of 
critical exponents and actual MC data. As a result, we have found that our 
values (77 = 1/10, 7 = 19/14, u> = 3/5) are consistent, whereas those of the 
RG theory (77 ~ 0.0355, 7 ~ 1.3895, uj ~ 0.782) are rather inconsistent with 
the MC data. 

4. Recent Monte Carlo data for complex zeroth of the partition function in 3D 
Ising model have been discussed (Sec. ^). The actual MC data suggest that, 
like in 2D Ising model, the critical behavior of the real part differs from that of 
the imaginary part. It can be explained reasonably by our exponents v = 2/3 
and uj = 1/2, but not by those (y ~ 0.6304 and ui ~ 0.799) of the conventional 
RG theory. Our best estimate of the critical exponent v from the MC data, 
i. e. 1/v ~ 1.518 or v ~ 0.659, is in a good agreement with the theoretical 
(exact) value 2/3. 
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